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Abstract 



The anisotropic Manev problem, which hes at the intersection of classical, 
quantum, and relativity physics, describes the motion of two point masses 
in an anisotropic space under the influence of a Newtonian force-law with 
a relativistic correction term. Using an extension of the Poincare-Melnikov 
method, we first prove that for weak anisotropy, chaos shows up on the 
zero-energy manifold. Then we put into the evidence a class of isolated 
periodic orbits and show that the system is nonintegrable. Finally, using the 
geodesic deviation approach, we prove the existence of a large non-chaotic 
set of uniformly bounded and coUisionless solutions. 
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Key words: Manev problem, anisotropy, chaos, geodesic deviation 



I. Introduction 

The Manev problem studies the dynamics of two point masses in a flat space 
for an attraction law given by a potential of the form ^ + ^, where r is the 
distance between particles and A,B>0 are constants. The correction term 
to the Newtonian potential provides a classical approximation to general 
relativity The anisotropic Manev problem, proposed by the first 

author in the mid 1990s [3], replaces the fiat space with an anisotropic one. 

The main reason for considering this question was that of seeking sim- 
ilarities between classical mechanics, quantum theory, and general relativ- 
ity. Its study was inspired by the anisotropic Kepler problem introduced by 
Gutzwiller in the early 1970s [Q. Gutzwiller aimed to find connections 
between classical and quantum mechanics. His interest in the anisotropic 
Kepler problem was aroused by an old quantum mechanical question related 
to a paper of Einstein [^]. 

Even though for a special class of (integrable) systems the Born-Sommer- 
feld-Einstein conditions seemed appropriate for describing a Coulomb 
limit of quantum theory, it was unclear how to find a classical approxima- 
tion for ergodic systems. Two qualities made the anisotropic Kepler problem 
suitable for addressing this issue: its chaotic character and its suitability to 
model various physical phenomena, as for example those encountered in cer- 
tain semiconductors [Q. Using a classical evaluation of the quantum Green's 
formula, Gutzwiller ^ found an approximate spectrum, in good agreement 
with previous quantum calculations. All these properties of the anisotropic 
Kepler problem raised hopes for the study of the Manev analogue, which also 
brought general relativity into the game. 

Some of the qualitative features of the anisotropic Manev problem prob- 
lem have already been studied. In a previous paper we proved the existence 
of an open, connected invariant manifold of uniformly bounded, collision- 
less orbits that lie within the negative-energy manifold [^. In 0] it was 
shown that the differential equations describing the anisotropic Manev prob- 
lem exhibit properties characteristic to all three aforementioned branches of 
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physics. That paper also brought arguments favoring the nonintegrabihty 
of the system by putting into the evidence a network of heterochnic orbits 
within the zero-energy manifold. Up to now, however, the existence of chaos 
and the nonintegrabihty of the system have not been rigorously proved. 

The first result of this paper shows the presence of chaos for a small set 
of solutions within the zero-energy manifold (see Section III). We will base 
our proof on an extension of the Poincare-Melnikov method developed in an 
earlier paper This method was used there to show the presence of chaos 
in several problems, including the restricted circular 3-body problem and the 
Gylden problem. The latter had been earlier identified as a degenerate case 
for applying the Poincare-Melnikov method [|lOl . In the 3-body problem the 
perturbation function, W = W{r, 9, t), where (r, 9) are the polar coordinates, 
appeared in its most general form. In the Gylden problem this function was 
isotropic, W = W{r,t). From the methodological point of view, we will 
now complete the applicative picture by offering a time-independent but 
anisotropic perturbation function of the form W = W{r,9), which is specific 
to the anisotropic Manev problem (see Section II). 

It is interesting to remark that different perturbations break up different 
symmetries, which affect the Melnikov integrals. In the general case both 
the homogeneity of time and the rotational symmetry are destroyed. This 
leaves us with two independent Melnikov conditions. In the Gylden problem 
the rotational symmetry is preserved (the angular momentum is constant), 
which leads us to one Melnikov condition. In the present case only the time 
homogeneity is preserved (through the conservation of the Hamiltonian) , so 
the only object of study is the other Melnikov condition (see Section III). 

Our second result proves that the anisotropic Manev problem has no other 
analytic integrals beyond the Hamiltonian function. To show this, we use a 



method designed by Poincare [11|, |12|, which we present in Section IV. The 



idea is to prove the existence of at least one isolated periodic orbit, which 
obstructs the occurrence of other analytic integrals. To apply Poincare's 
method we need to write the perturbation as a Fourier series. This leads to 
some technical difficulties, which we overcome in Section V. In Section VI 
we prove the existence of many isolated periodic orbits in the hypothesis 
of weak anisotropy. We then show that this implies nonintegrabihty. The 
existence of periodic orbits, however, is an interesting result in itself, which 
raises hopes for applying Gutzwiller's formula to this problem in the context 
of semiclassical studies of the atom. 
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In the last part of the paper we are concerned with the quahtative be- 
havior of colhsionless orbits for which the two particles stay close to each 
other. In Section VII we present the geodesic deviation method [|T3|, which 



provides a criterion for excluding the presence of chaos. The idea is to trans- 
form the equations of motion into some equivalent ones. The new equations 
put the solutions of the given system in a one-to-one correspondence with 
the geodesies of a certain Riemannian manifold whose metric is determined 
by the initial equations of motion. Measuring the deviation of the geodesies, 
we can draw conclusions about the behavior of the initial solutions. If the 
geodesies diverge, chaos may appear; if they converge, chaos is ruled out. In 
Section VIII we show that there exists a positive-measure invariant mani- 
fold of orbits that fail to encounter collisions but for which the particles stay 
close to each other. Using the geodesic deviation method we then prove that 
the corresponding geodesies converge, thus showing that the orbits of this 
manifold are not chaotic. 



II. Equations of motion 

The (planar) anisotropic Manev problem is described by the Hamiltonian 

^(q,P) = ^P - / 2 2 ~ nr?^c? ' 



where 6 > and /i > 1 are constants, q = {qi,q2) is the position of one body 
with respect to the other considered fixed at the origin of the coordinate 
system, and p = {pi,P2) = q is the momentum of the moving particle. The 
constant n measures the strength of the anisotropy; the larger n, the higher 
the difference between the weakest and the strongest directions of attraction. 
For /i = 1 we recover the classical (isotropic) Manev problem of which we 
have a complete quahtative understanding [0. Therefore, if fi is only 
slightly larger than 1, we are dealing with a weakly anisotropic case. Unlike 
in the isotropic problem, the angular momentum K{t) = p{t) x q(t) of the 
anisotropic case is not a first integral of the system. This is because the 
rotational invariance of the Hamiltonian breaks for /i > 1. Therefore we 
expect to encounter richer dynamics in the anisotropic case. 

In the first part of this paper we will be interested in aspects related to 
weak anisotropy, i.e., in values of ^ > 1 that are close to 1. To put into 
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the evidence the perturbative character of the anisotropy with respect to the 



isotropic problem, we introduce the notations e = fi — 1 and r = ygf + gl- 
For n > 1 close to 1, the quantity e > is small, so the Hamiltonian (1) can 
be viewed as describing a perturbation of the Manev problem by writing 



n{q, P) = lp'-l-^ + e(^ + ^] cos^ 0^no + e>V(r, 6). (2) 



As in ||10[, we now consider the parabolic solutions of the unperturbed 
problem (defined by the Hamiltonian Tio) that lie on the zero-energy mani- 
fold. These orbits, which play the role of homoclinic solutions corresponding 
to the critical point at infinity {r = oo, f = 0), satisfy the equations 



where /c 7^ is the constant angular momentum of the unperturbed problem; 
the negative (positive) sign corresponds to t < (t > 0). From (3) we get 
that 



±t = !s^^f+L^2r-{k^ -2b) + constant 
9 = ±2 arctan "^^^^^.2^26 ~'~ constant. 



(4) 



Let 

R = R{t) and 6 = e(t) (5) 

be the expressions of r and 6 with respect to t, obtained by inverting equa- 
tions (4) and assuming that -R(O) = r^m = k'^/2 and G(0) = vr. We will 
not need their expressions; as in 0, we only retain the information that R 
is even and that is odd. 



III. The Melnikov integral 

We now consider the problem defined by the Hamiltonian Ti given in (2). 
We will call homoclinic manifold the set of solutions of the unperturbed 
system that are doubly asymptotic to the point r = 00, r = 0. For k ^ 0, 
the two-dimensional homoclinic manifold is described by r = R{t — to) and 
= Q{t — to) + 60, with arbitrary constants to and ^^o- 
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It is interesting to remark that in this problem we are deahng with nega- 
tively and positively asymptotic sets to the critical point r = oo, f = and 
that this point would correspond to the "critical point at infinity," obtained 
with the help of McGehee's transformation |15| 



dt = —ds. (6) 



3j 3j 



But instead of using McGehee's technique, we will apply here the method 
introduced in It was shown there that under natural assumptions imposed 
on the perturbation function }V{r,6), one can guarantee the existence of 
smooth solutions approaching the point r = oo, r = Ofort^ ±oo. It was 
also proved that the Poincare map leads to infinitely many intersections of 
the stable and unstable manifold for the corresponding fixed point. 

In our case the perturbation arising from a weak anisotropy vanishes as 
r —>■ oo. This happens because yV{r,6) ~ 1/r, which is exactly condition 
(18) in [0. So we can write the Melnikov condition as in [9], with the only 
amendment of dropping the time dependence. Therefore the two Melnikov 
integrals become 



r+oo 

MM 



and 



-^^^^^ dwm),m+do) ^ Q^^^9w(i?(t), 0(t) + ^o) 



dr 09 



dt = 
(7) 



M.(.o)^/:;™)^^^^0. (8) 
Since W vanishes as t ^ ±oo, the first Melnikov condition takes the form 

M,(«„)=r5M^(%^W±M,«.o. (9) 

J — oo ut 

Ml is identically zero since the perturbation function W is independent of 
time. This simplifies our discussion because, unlike in the general case, we 
only need to find the zeroes of the integral in (8). 

It is significant to remark, and easy to check, that the previous conditions 
can also be written in terms of the first integrals of the unperturbed problem 
as 

/ + 00 
{no,w}{...) dt = o (10) 
-oo 
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and 

/ + 00 
{K,W}{...) dt = 0. (11) 
-oo 

Notice that Mi{9q) is identically zero. This resembles a results obtained for 
the Gylden problem p, |10[ and is related to the symmetries of the problem. 
In the Gylden problem the perturbation is independent of the angle 9 but 
depends on time. This means that the perturbation alters the homogeneity 
of time, so the Hamiltonian is not an integral of motion anymore, but it 
leaves the rotational invariance intact; thus the angular momentum is still 
conserved. Therefore there is only one condition here, given by relation (9). 
On the other hand in the anisotropic Manev problem the anisotropy of the 
space alters only the rotational symmetry but not the homogeneity of time. 
This also leads to only one condition, given by relation (8). 

For the anisotropic Manev problem the Melnikov condition M2 takes the 
form 

M2(^o) = sin[2(e(t) + ^0)] [^-^^ + -^^jdi = 0- (12) 
Using some trigonometry the integral can be written as 

M2(^o) = h cos 260 + h sin 2^0, (13) 
where /i and I2 are defined by the relations 

^i = /-c^(i^P)+w)sm2e(t)rft 

+00 / 1 I b 
-00 I R{t) "t" 



/-r(7^ + 7^)cOs2e(t)rft. 



(14) 



Since R and are even and odd functions of time, respectively, the integrand 
of Ji is an odd function. Therefore Ji = 0, and finally M2 can be rewritten 

as 

M2(0o) =/2sin2^o. (15) 

Since we found simple zeroes of the Melnikov function, the following result 
holds. 

Theorem 1 In some invariant set contained in the zero-energy manifold, 
the weakly anisotropic Manev problem exhibits chaotic dynamics. 
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This type of chaotic behavior, which usually takes place for a small set 
of orbits, is induced by the infinitely many intersections of the stable and 
unstable manifolds for the Poincare map associated to the critical point at 
infinity. Let us notice that the Poincare-Birkhoff-Smale theorem does 
not directly apply to this situation, which is degenerate. However, the exis- 
tence of Smale horseshoes and of positive topological entropy is not restricted 
to hyperbolic equilibria. This phenomenon is also encountered in nonhyper- 
bolic for example those arising when dealing with area-preserving 

perturbations fl^. Thus, Theorem 1 adds to the class of results describing 
chaotic dynamics near degenerate equilibria. 



IV. Poincare's method 



In this section we will present a classical result of Poincare |Tl| , |]T2| in con- 
nection with the nonexistence of additional analytic integrals in Hamiltonian 
systems. Based on the investigation of long periodic solutions, this criterion 
is suitable for proving the nonintegrability of the anisotropic Manev problem 
for weak anisotropy, as we will see in Section V. 

First we need to write the unperturbed system, i.e. the (isotropic) Manev 
problem, in terms of action-angle variables |T^. The action variables are 
given by 

[K = qiP2 - q2Pi, 

where h is the energy constant and K is the angular momentum. These vari- 
ables are defined for h < and > 2b. The associated angular frequencies 
are 

^ (/+VX2-2fe)3 

, , ^ K 

With the help of the action-angle variables, the unperturbed Hamiltonian 
Tio defined by relation (2) can then be written as 

Ho{I,K) ^ 



2{I + VK^-2by' 
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so the unperturbed equations of motion take the form 



= ^ = 



(16) 



Recall that the eigenvalues A of the monodromy operator of a T-periodic 
solution are called characteristic multipliers and that the numbers a defined 
by A = exp{aT) are called characteristic exponents ||16|, |18[. We would 
like to show the nonexistence of other analytic integrals, independent of the 
Hamiltonian Hq. Our main theoretical tool in proving the nonintegrability 



of the perturbed system is the following result, proved in |TT 



Theorem 2 Suppose that the two degree of freedom Hamiltonian system with 
Hamiltonian function H has two first integrals H and F that are independent 
along a periodic solution. Then four characteristic exponents vanish. 

One exponent vanishes because the system is autonomus. The second 
vanishes due to the presence of the first integral H. If the remaining expo- 
nents are different from zero, the periodic solution is called nondegenerate. 
It is well known that nondegenerate periodic orbits are isolated |1T2| . 

To establish the existence of a large number of isolated periodic orbits, 
suppose that for / = /° and K = K^, the frequencies ui and uk of the 
unperturbed problem are commensurable and that ui ^ 0. Then the per- 
turbing function W{I, K,uJit,ujKt + A), defined by W given in (2) and then 
transformed with the help of the action-angle variables, is periodic in t and 
has period T. Consider its average 

Tr(/°, \) = ^ r W{I, K, ujt, UKt + X)dt. 
I Jo 



Then the following theorem, also due to Poincare ||ri|, establishes the exis- 
tence of isolated periodic orbits. 

Theorem 3 Assume that the following conditions are satisfied: 

(1) The Hamiltonian is nondegenerate at the point I = I^,K = , 

(2) for some X = X* the derivative dW/dX = 0, but d'^W/dX'^ ^ 0. 
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Then, for small e ^ 0, the perturbed Hamiltonian system has a T-periodic 
solution that depends analytically on the parameter e and for e = coincides 
with the periodic solution 

I = I\ K = K°, (j) = uJiT, e = ujKt + X* 

of the unperturbed system. The two charateristic exponents ±a of this solu- 
tion admit a convergent series expansion in power series of ^/e, 

a = aiy/e + + a^ey/e + . . . , (17) 

where 

problem. Using Theorem 2, one can show the dependence of the functions 
i^o and Fq on the set of unperturbed tori I = Iq, K = Kq that satisfy the 
conditions of Theorem 3 and the relation 

So the following result holds. 

Theorem 4 Under the hypotheses of Theorem 3 and relation ( 19), the Hamil- 
tonian system (16) does not admit a first integral F independent of H that 
can be written as a formal power series of the form I]s>o Fs{I, K, (p, 6)e^ with 
analytic coefficients Fg. 



V. Fourier series of the perturbation 

In order to apply these results to our perturbed problem, we first need to 
write the perturbation function as a Fourier series. For this let us first notice 
that r can be expressed in parametric form with respect to time. Following 
the same procedure used in [|l^] for the Kepler problem, we can write 

/ 1 r rdr 

^~ \2\h\J l-r'^ + J- _ K'-2b- 
' ' ' V I'*! 2|/i| 
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With the substitutions 

A = ^— and A^e^ = ( ^ri] - ^ , 
2\h\ V2|^l/ 2|/i| ' 

the above relation takes the form 

rdr 



^2g2 _ _ j\y 



With the change of variables r — A = — Aecosr/, the integral assumes the 
simpler expression 

t = A 2 y (1 — e cos vi)dri = A 2 {rj — e sin rj) + constant. 

If we choose the time in such a way that the addition constant is zero, then 

3 

r = A(l — e cos rj) and t = A^ (^rj — e sin rj). 
We can now write the constants A and e in terms of the action variables as 

2\h\ ^ ' (/ + _ 26)2 

Recall that the Fourier series expansion of a function F has the form 

F{ri) = V Aie'^", where A = — F{ri)e-'^'^dT], 
i=-oo 27r Jo 

and that that the Bessel function Jm{z) is given by 



Itx Jo 

It is now easy to show that 

A 
r 



1 

— / cos(mu — z sin u)du = (—l)'^J-m(z). 
In Jo 



1 + 2 ^ J; (/e) cos Ir], 
1=1 



where rj = ujjt. Therefore 



CO 00 

— = 1 + 4^ Ji(/e) cos/?7 + 4 ^ J i{le). J m{'me) cos Irj cos mrj. 



r 



1=1 Lm=l 
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Using some trigonometry we obtain that 

^ = 1 + 4 Ei=i Mle) cos lrj{l + m)+ 
2Ei^=i Ji{le)Jm{me)[cosr]{l + m) + cosr]{l - m)]. 
More computations lead us to the relation 

— = 1 + Y.iMle)f + iMle) + y + a) cosrjl - 2A, cosrj, 



1=1 



where 



Ai^ ■Ja{ae)J(3{pe) and Bi ^ Yl Ja{ae)Jp{pe). 



a,l3 
a+0=l 



We can write the the perturbation function as a series of the form 

oo 

W(I, K, ujit, uxt) — + ^ Q cos lujit}j cos^ ojxt^ 

where 



1=1 



C, 



Ah 



Ji{le){^ + l) + ^ + Bi\ for 1^1, 



Ai 
2 



and 



Ci = %lMle)(j + l) + Y + Bi\-2A„ 



h °° 



1=1 



Using some trigonometry, the above relation becomes 

W{I,K,uit,UKt) ^ D + - Y Ci^rnCos{lu;i + 2mu}K)t, 

^ (=0 

m=±l,0 



where C/,^ = Ci for all l,m ^ 0, Co,i = D, and Cq-i = 0. 
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VI. Periodic orbits and nonintegrability 



We can now apply the method of the previous section to the anisotropic 
Manev problem. As we are going to show, all the assumptions of Theorem 3 
are satisfied, so for weak anisotropy the existence of isolated periodic orbits 
follows. This leads to the following result. 

Theorem 5 For small e ^ 0, the weakly anisotropic Manev problem has a 
T-periodic solution that depends analytically on the parameter e. For e = 
this solution coincides with the periodic solution 

I^I^, K=K\ 4)^0JiT, e^UKt + X* 

of the unperturbed system, where A* can take the values 0, |, tt, Moreover, 
the two charateristic exponents ±a of the solution admit a convergent series 
expansion in ^/e, given by (17), where is defined by (18). 

Proof. First we will show that the assumptions of Theorem 3 are satisfied. 
It is easy to see that the unperturbed Hamiltonian is non-degenerate. Indeed, 




dIdK 



dKdl dK^ / ^ y V I V I 

We are left with computing the average W of the perturbation function W for 
values of the action variables whose frequencies of the unperturbed problem 
are commensurable. Let ujk — ^^i- Then we can write the perturbation as 



Ci^rn cos 

^ 1=0 

m=±l,0 



I + 2m- ] ojit + 2mA 



It is clear that the only A-dependent terms that survive after averaging are 
the resonant ones, i.e., the terms with the property that I + 2m^ = 0. It 
is easy to sec that we have to consider only the terms with m = ±1. This 
implies g = ±1, ±2 and then / = ±2p, ±p. These terms are of the form 

C±2p,±i cos (2A) and C±p,±i cos (2A) , 

so 

dW d'^W 

G sin 2 A and — — =2Gcos2A, 



dX dX^ 
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where 

G = 2 (C_2p,i + C2p,i + C-p^i + Cp^i) . 

This means that for A* = 0, |,7r, the second condition of Theorem 3 is 
satisfied if G is not identically zero. 

We will now show that, in general, G cannot be identically zero. Indeed, 
notice first that G depends on A and is analytic in A. But G cannot be 
identically zero beyond a discrete set of values of A unless the Bessel function 
is itself identically zero. However, the Bessel function has only a discrete 
set of zeroes. Therefore, generically, G cannot be identically zero, so the 
existence of isolated periodic orbits follows. This completes the proof. 

For the anisotropic Manev problem, we are now able to prove the nonex- 
istence of first integrals that are independent of H and analytic in the pa- 
rameter e. In more formal terms, the following result holds. 

Theorem 6 For small values of e, the weakly anisotropic Manev problem 
does not admit a formal first integral F independent of H that can he writ- 
ten as a power series Y.sXd Fs{I, K, cf), 6)e^ whose coefficients Fg are analytic 
functions. 

Proof. According to Theorem 4, the only thing we need to prove is that 
inequahty (19) is satisfied. Then using Theorem 5 the result follows. In our 
case inequality (19) can be written as 

If we substitute = |w/ (where cjj ^ and p 7^ 0) we obtain that (19) is 
always verified except for the solutions of the equation 

But this equation has at most five distinct solutions, which do not affect the 
outcome since, if eliminating them, the remainig set is still dense in the real 
fine. This completes the proof. 
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VII. Geodesic deviation 



In this and the following section we will show that if the particles do not col- 
lide but stay close enough together, the motion is not chaotic. This is some- 
what surprising since one expects that, as in celestial mechanics or atomic 
physics, chaos appears because of near-collision approaches. Of course, our 
result does not exclude the possibility that chaos is determined by only brief 
passages close to collision. 

Our proof is based on a local criterion due to Szcz§sny and Dobrowolski 
The idea is to rewrite the equations of motion in terms of the geodesic 
equation with respect to the Jacobi metric and to measure the local tendency 
of the geodesies to diverge or converge. If they diverge, chaos may appear, 
but if they converge, chaos is impossible. 

Let us first outline the ideas that lead to the above mentioned criterion. 
In general, for a Hamiltonian function 

where g^^ are the components of the covariant tensor corresponding to the 
Riemannian metric Qijdq^ ® q-', with pj = Qijc^^, the equations of motion can 
be written as the Euler-Lagrange system 

ci' + T\^ci>^ci' = -g^'dkUiq), (20) 

in which F'^^ is the Levi-Civita connection for the metric Qij. The energy 
relation ?^(q, p) = h foliates the phase space in codimension-one manifolds. 
Let us fix an h and define an open set D/j of the configuration space as 
D/i = {q| U{q) < h}. It can be shown that the orbits of (20) are orthogonal 
to the boundary of D/j. Let us define the Jacobi metric g^^ = 2{h + U)g^^ 
and the parameter s such ^ = 2{h + U). It can be shown that on the set 
Dh, equations (20) are equivalent to the geodesic equations 

ds^ + ^ ds ds ~ 

where r*^-^ is the Christoffel symbol with respect to the metric g'^K Note that 
5'ij(q(s))^^ = 1. Equations (21) can be further written as 

Vuu = 0, 
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where u is the tangent vector to the geodesic and V is the covariant deriva- 
tive with respect to the Jacobi metric. Aiming to understand the relative 
motion of geodesies in a domain inside D/^, we take a curve transversal to 
the geodesies, then take a copy of it by transporting each point of the initial 
curve along the geodesic, and measure in this way the local relative deviation 
of the geodesies. It can be proved that, in classical form, the measure of the 
deviation is given by the solution of the Jacobi equation 

^ - ^>C(s)., (22) 

where x is the variable that measures the deviation and /C is the Gauss 
curvature at x{s). If /C > the geodesies approach each other, whereas if 
/C < they diverge. For 2-degree of freedom systems the curvature is given 



by the formula |]13 



_ g^^{d,V)id,V) + ih-U)AU 
^ - 4(/. - f/)3 ' ^^^^ 



where A = df + 82 is the Laplacian operator. In [13| the following result is 
proved. 

Theorem 7 IfJC given by formula (23) is negative, the geodesies diverge and 
chaos may appear; if K, is positive, the geodesies converge and chaos cannot 
take place. 

We will further use this result to investigate the possibility of encountering 
chaotic motion in the neighborhood of collisions. 



VIII. Absence of chaos near collision 

In order to apply Theorem 1 to the equations of motion of the anisotropic 
Manev problem, let us first show that the hypothesis that leads to the 
geodesic deviation equation (22) is satisfied. For this we need to make sure 
that the geodesies do not reach the boundary of the set D/i, in other words 
we must isolate a set of solutions with this property. This will be achieved 
while proving the following result, which shows that collisionless solutions 
that keep the particles close together cannot be chaotic. 
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Theorem 8 There exists a positive-measure set of orbits that are free of col- 
lisions and for which the distance between particles does not exceed a certain 
uniform bound. Within this set, chaotic motion is ruled out. 

Proof. In the case of the anisotropic Manev problem, the expression of the 
kinetic energy shows that the metric is fiat and that gij — 1 ior i — j and 
Qij — ior i j. Therefore the curvature JC given by formula (23) takes the 
form 

_ \VU\' + jh - U)AU 
^ 4{h^^ ' (2^) 

where U(q) = , ^ j^-t- Let us now see that there exists a set of 

solutions, whose projection in the configuration space lies inside D/j. These 
will form an open set of uniformly bounded and collisionless solutions. 

For this consider the analytic diffeomorphic transformation of the depen- 
dent variables, 

r = (ql + qi)^ 
ip = arctan(g2/q'i) 
v^r = qipi + q2P2 
k u = r^(fi = qip2 - ?2Pi, 
and the analytic diffeomorphic transformation of the independent variable 

dr = r-^dt. (25) 

Notice that r and (/? are polar coordinates, whereas v and u represent the 
rescaled radial and tangential components of the velocity. In these new vari- 
ables, the equations of motion take the form 



r' — V 



(26) 



v' — 2hr + {n cos^ ip + sin^ ^)~^^^ 
(fi' — r~^u 

_ u' = ^r~^[(/iCOS^ cp + sin^ </?)~^/^ + 26(//cos^ (p + sin^ <p)~'^] sm2(p 
and the energy relation ?i(q, p) — h becomes 

u"^ + v'^ - 2r{ii cos^ (p + sin^ p})'^^"^ - 2b{id cos^ (p + sin^(p)~^ = 2hr'^. (27) 

The new dependent variables (r, v, (p, u) G (0, oo) x IR x 5"^ x H are functions 
of the fictitious time r, so the prime denotes from now on differentiation with 
respect to this new independent variable. 
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Let us first notice that there exist two constants, mi and 1712 with < 
mi < m2 < 00, such that along any solution (r, f,^?,^) of system (26), the 
function of r given by // cos^ + sin^ satisfies the relations 

mi < jjL cos^ Lp + sin^ Lp < m2, 

for all r for which ip is defined. Consequently there exist two constants. Mi 
and M2, with < Mi < M2 < 00, such that 

Ml < {11 cos^ cp + sin^ cp)'^^^ < M2 (28) 

for all r. From the first two equations of system (26), we obtain the nonho- 
mogeneous second-order equation 

r" - 2hr = (// cos^ cp + sin^ ¥^)"^/^ (29) 

which represents a forced harmonic oscillator. We will further ignore the last 
two equations of system (26) and retain only the information that the forcing 
function (/icos^v? + sin^v?)"^/^ is uniformly bounded from above and away 
from zero from below, as relations (28) indicate. 

We will further restrict our analysis to the invariant manifold of solutions 
of negative energy, h < 0, which obviously exists according to relation (27). 
Solving the homogeneous equation r"—2rh — and then applying to equation 
(29) the method of variation of parameters, we obtain the general solution 
of equation (29) in the form 

^(^) = (^1 + 7^ /o -(^-)^ da) sin(y=2^r) 
+ (^2 - ^ lo -yv^-)^ da) cos(v^r). 

\ V ^1 n COS-' (p(a)+sm'' (p(a) ' 

If, for every solution of constants Ci and C2, we apply to each of the above 

integrals the intermediate- value theorem and then use some trigonometry we 

obtain ^ 

r(T) = (cf -I- c|)2 cos(a;T — a;o) — 

sin2(\/-2/iT) cosp-{\J -2hT) 

(30) 

2h^/lCOs2 ly5(Tl)+sin^ (p( taU\) 2/l^ ^COS^ ip{T-l)-\-S\T? •f(T2) ' 

where ti and T2 belong to the interval (0, r), = \^—2h and coq = arctan — . 
Relations (28) and (30) allow us to draw the conclusion that for any solution 
{r,v,ip,u) of system (26), the inequalities 

^cf + ci cos(^r - ^0) - ^ < r(r) < ^c? + 4 cos(a;r - cuq) - f (^1) 



17 



take place for all r for which the solution is defined. 

Let us now fix some /iq < 0, as close to as we like, and define the set 
A(/io) of solutions formed by the union of all sets A^^, for h in the interval 
(— oo,/io); where each set Ah contains all the solutions given by (30) that 
satisfy the inequality 



Ml 



—1— f- <r 



Obviously, this is an open, nonempty, and connected invariant manifold of 
system (26). Prom (31) it follows that, for every solution of this manifold, r 
is positive and bounded, therefore the orbits are coUisionless and bounded. 
Let us denote by A(h^ the projection of A(/io) into the configuration space. 
For fixed, suitable choices of ci,C2, and h provide solutions with qy and 
52 small enough such that V{c\) < ho. Moreover these choices can be made 
such that r is as close to as we like but still positive. Obviously, the set 
of the corresponding solutions has positive measure. Let us denote by T{ho) 
the corresponding projection of the set of such solutions in the configuration 
space. Obviously T{ho) is contained in A(/io) as well as in D/j^. This shows 
that the hypothesis leading to the criterion of the previous section is satisfied. 
A straightforward computation yields 

^ _ 3{h-U){^^q^+qi) _ [l-8b(fc-t/)](/x^g?+gi) , 4bHp,^q^+q^) 4b{,,^q^+ql) 

The first term is always negative and the second can also be negative, de- 
pending on the fixed value of h. But since the last two terms are positive and 
the powers at their denominators exceed those of the denominators of the 
first two terms, for sufficiently small q\ and ^2, the curvature /C is positive. 
We can always define r(/io) such that the configurations corresponding to 
the solutions with corresponding and qi are contained in r(/io)- There- 
fore, according to Theorem 7, the set of coUisionless solutions for which the 
bodies stay close enough together, does not exhibit chaotic behavior. This 
completes the proof. 
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